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Abstract 

The form factor of baryons is estimated using 3.46 ■ 10^ hadronic Z decays 
collected by the DELPHI experiment between 1992 and 1995. Charmed 
baryons fully reconstructed in the pK~7r~^, pKg, and A7r"'"7r"'"7r~ modes, are 
associated to a lepton with opposite charge in order to select A° K^l'Vi 
decays. From a combined likelihood and event rate fit to the distribution of the 
Isgur-Wise variable and using the Heavy Quark Effective Theory (HQET), 
the slope of the 6-baryon form factor is measured to be : 

p2 = 2.03 ± 0.46 (stat)^°:^g (gyg^)_ 

The exclusive semileptonic branching fraction i?r(A° A^l^ui) can be de- 
rived from and is found to be (5.0lQ;g (stat)^}'2 (syst))%. Limits on other 
branching fractions are also obtained. 
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1 Introduction 

The knowledge of the S-meson form factor has recently improved thanks to a wealth 
of new experimental results, as reported for instance in references [1-6]. Semileptonic 
decays of B mesons into D and D* final states can be understood in the context of the 
Heavy Quark Effective Theory (HQET), where the four form factors remaining when the 
lepton mass is neglected can be expressed in terms of a single Isgur-Wise function ^, 
which will be defined in section 2. 

In this paper, the semileptonic 6-baryon decay A° A+/~z/; (with /~ = e~ or is 
investigated ^, where the is fully reconstructed from its decay modes into pK~7r~^, 
Att+tt+tt", and pK^. The heavy quark symmetry relates form factors of the transition 
A° — > A+Z~P; to a new single Isgur-Wise function ^s, as explained in [7-9] and predicts 
its absolute value when the final A+ is at rest in the Aj, frame. After a summary of 
the heavy quark formalism in section 2, and a description of the relevant parts of the 
DELPHI detector in section 3, the selection of the different A+-lepton channels candidates 
is described in section 4. The dominant contribution to the A^l'Pi final state is expected 
to come from the A^ baryon, and the relevant contaminations from S-mesons, from other 
6-baryons, or from other hadronic final states with additional pions are investigated in 
sections 5 and 6. In section 7 a direct determination of the semileptonic branching fraction 
is presented, while the fit to the distribution of the Isgur-Wise variable w is described in 
section 8. The semileptonic branching fraction and the w distribution are then combined 
into a single measurement of the slope parameter of the 6-baryon form factor, assuming 
the validity of HQET predictions. 



2 The semileptonic decay form factor of b-baryons 

A complete description of the form factor formalism in semileptonic decays and the- 
oretical predictions can be found in [10-12]. The form factors are functions of the four- 
momentum transfer in the transition, with = {pi +pp;)^ where pi and pp^ are the 
charged lepton and neutrino four-momenta, respectively. Isgur and Wise introduce the 
dimensionless variable w, scalar product of the four- velocities of the A° and A+: 

w = Va^- Va^ = {ml^ + - q^)/ (2mA^mAj. (1) 

The hadronic current in the weak decay of a beauty baryon (J^ = 1/2+) to a charmed 
baryon [J^ = 1/2+), as in the transition A^ A'^W^, involves six form factors, three 
vectors and three axial-vectors Gj. In the decay A° — > A+Z~P;, the variable w ranges 
from 1 (highest transfer, final A+ at rest) to a value close to 1.44 (smallest transfer, 
= mf). Among the six form factors which can contribute to the semileptonic decay of 
= 1/2+ baryons, only Fi and Gi survive in the limit of infinite b and c quark masses, 
the HQET limit. In addition, they are equal and can be expressed in terms of a single 
function ^b(w). 

The differential decay width of the transition A° — > A'^l~ui can be obtained from [13] 
in the approximation where ninght/mQ terms are neglected, where mught is the mass of 
the light quark system and uiq stands for the heavy quark: 

^^GK{w)eB{n^), (2) 

^The subscript M has been added to indicate that this form factor applies to meson decay, while the subscript B will 
be used for baryons. 

^The notations AjJ and AjJ" will implicitly stand for both the baryon and antibaryon, with the proper inversion of the 
signs of the lepton and of the A^ decay products. 
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where the constant G is: 

G = \^ch? m\j^ with r = mxJmK^, (3) 

and the kinematical factor K{w) is: 

K{w) = P [3w{l - 2rw + r^) + 2r{w^ - 1)] with P = mA^Vw^ - 1. (4) 

The Isgur-Wisc function C,b{w) will be studied in the present paper. 

In i?-meson decays, another function ^Miw) describes the semileptonic transitions 
B — > {D, D*)l~Pi, and a Taylor expansion is usually assumed for this function: 

^mH = eM(l) [1 - pliw - 1) + 0{{w - If)] , (5) 

with ^m(I) = 1 in the HQET limit. The quadratic terms are constrained by dispersion 
relations [14-16]. Taking into account {mught/Tnc) corrections and perturbative QCD 
effects [10-12], the value of Cm(1) is modified in the 5-meson decay channel to {m(1) = 
0.91 ± 0.04. Several experimental determinations of the mesonic form factor ^m{w) have 
been performed in the channel B'^ — > D*'^l~Vi. Different fits are performed in [4] which 
quotes for a constrained quadratic fit: 

= 1.22 ±0.14 (stat). (6) 

Another recent determination is given in [5], and earlier measurements of p\j were 
performed by [1-3]. 

Many different parametrisations have been proposed for the baryonic function ^b{w), 
as given in references [13,17,18]. The simplest one is chosen which remains positive in 
the physical w range: 

eBH=eB(l)xexp[-p2(^^-l)]. (7) 

The linear and quadratic coefficients in the Taylor expansion of this function $,{w) are 
(within errors) in the domain allowed by the dispersion relation constraints evaluated 

in [14-16]. 

The flavour independence of QCD implies that, as in the usual isospin symmetry, 
isiX) = 1- It is shown in reference [13] that the corrections in mught/mQ to this result 
vanish at first order for baryons (as expected from the general result of [19]), and remain 
small at higher order, while this is not true for mesons. The relations between the six 
baryonic form factors and the Isgur-Wise function ^ are however slightly modified by 
the small mught/mQ corrections evaluated in [13]. The perturbative QCD corrections 
are smaller than for mesons and will be neglected, as explained in [20]. In sections 7 
and 8 of this paper where the observed exclusive semileptonic branching fraction is used 
to infer p^, the finite mass corrections, as given by [13], are included into the relation 
between the Isgur-Wise function ^b(u') and the semileptonic width (equation 2). These 
corrections are evaluated with a h quark mass = 4.844 GeV/c^, and a charmed quark 
mass rric — 1.35 GeV/c^, the numerical values chosen in [13]. 

3 The DELPHI detector and the simulation 

The DELPHI detector and its performance have been described in detail in [21,22]. 
In the barrel region, a set of cylindrical detectors, with the z coordinate axis oriented 
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along the 1.2 T magnetic field (and the beam) direction, allows the tracking of charged 
particles. The silicon Vertex Detector (VD), with an intrinsic resolution of 7.6 /xm in the 
plane transverse to the beam axis and 10-30 //m along the z axis, consists of three layers. 
The innermost and outer layers were replaced by double- sided silicon microstrips for the 
1994-1995 data taking period. The Inner Detector (ID) extends between radii of 12 cm 
and 28 cm and gives 24 spatial measurements. The Time Projection Chamber (TPC) 
provides up to 16 points between 30 cm and 122 cm. The Outer Detector (OD), at a 
radius of 197 cm to 206 cm, consists of 5 layers of drift cells. In the plane orthogonal to 
the beam direction, the extrapolation accuracy at the primary vertex of hadronic charged 
particles is found to be -^20^ -|- 65^/p^ A^m [23], where pt (in GeV/c) is the momentum 
of the particle transverse to the beam axis. 

The identification of electrons relies on the electromagnetic calorimeter in the barrel 
region (high density projection chamber HFC), with a relative energy resolution of 6.5% 
for electrons at a momentum of 45 GeV/c. Within the HFC acceptance, electrons with a 
momentum above 3 GeV/c are identified with an efficiency of 77%. The probability that 
a pion be misidentified as an electron is below 1%. The muon identification relies mainly 
on the muon chambers. The selection criteria used in this work ensure an efficiency of 
identification of 77% for a misidentification probability of 0.8%. 

The identification of protons and kaons relies on the DELFHI algorithms which take 
into account the information provided by the Ring Imaging Cherenkov (RICH), and the 
dE/dx in the Time Projection Chamber (TPC). The liquid and gas radiator signals of 
the RICH arc used when they were present, mostly in the 1994 and 1995 data samples. 
A neural network program is used for the identification of charged particles in the 1994 
data, while in the other data samples, a simple combination of the RICH and TPC 
measurements are considered. The proton and kaon tracks are required to be in the 
barrel region, with a polar angle to the z axis fulfilling |cos^| < 0.74. The efficiencies 
associated to proton and kaon identification have been obtained from the simulation, 
and corrected for the small differences between the simulation and the data. A check 
of the efficiency estimate was performed on dedicated samples of real data, namely A 
(for the proton) and reconstructed D*^ in the Kttt: channel (for the kaon): a good 
agreement was found with the simulation in the whole momentum range. The overall 
proton identification efficiencies are (24±4)% (1992 data), (21.5±4)% (1993 data), when 
the liquid RICH was not operating, and (42 ±2)% in 1994 and 1995. The fraction of 
pions misidentified as protons is approximately 5% above 3 GeV/c for good operating 
conditions of the RICH. 

Special samples of events for each potential source of A+-lepton final states were gener- 
ated using the JETSET 7.3 Parton Shower program [24], with a A^ lifetime set to 1.6 ps. 
The generated events were followed through the detailed detector simulation DELSIM. 
The background events without a true Ac were selected from the general DELPHI simu- 
lated sample of qq events at the Z. The same sample was used to estimate the fake lepton 
background. These events were then processed through the same analysis chain as the 
real data. A reweighting of the simulated events, which were generated with a constant 
value ^b(w) = 1 of the form factor, allows the observed distributions to be predicted for 
all variables, and the slope parameter to be tuned in order to reproduce the data. A 
reweighting was also applied to match the measured A^ lifetime of 1.23 ps [25]. 



4 



4 Event selection 

4.1 The sample of hadronic events 

Hadronic Z decays collected between 1992 and 1995 were used. The centre-of-mass 
energy was required to lie within 2 GeV of the Z mass. Charged particles were required 
to have a measured momentum between 0.1 GeV /c and 50 GeV/c, a relative error on 
momentum less than 100%, a track length larger than 30 cm, and a distance of closest 
approach to the interaction point smaller than 5 cm radially, and smaller than 10 cm 
along the beam axis. Neutral particles were required to have an energy between 1 GeV 
and 30 GeV, and a polar angle between 20 and 160 degrees. They were assigned a zero 
mass. 

Hadronic events were then selected using the previous set of charged particles with a 
momentum above 0.4 GeV/c. Five or more charged particles were required, carrying a 
total energy (assuming them to be pions) of more than 0.12-y/i. A total of 3.46 million 
events has been obtained. In the following sections, events will be selected from this 
sample which contain candidates for both a Ac and a lepton. 

4.2 Lepton selection and b-tagging 

The lepton candidates had to satisfy the appropriate identification criteria, and a 
method relying on a neural network was used for the electron identification. The charged 
leptons were required to have a momentum larger than 3 GeV/c and a transverse mo- 
mentum pt with respect to the A+ candidates (defined in sections 4.3 to 4.5 below) larger 
than 0.6 GeV/c. Additional criteria were then introduced to purify the A° sample: 

• The mass of the A+-lepton system was required to be larger than 3 GeV/c^ (for 
electrons), or 3.2 GeV/c^ (for muons). This selection reduces the potential contribu- 
tions from semileptonic 5-meson decays with baryons in the final state and from the 
A+X/~z/i final states of A° decays. The lower momentum selection for the electron 
preserves the efficiency, given a lower resolution in this channel. 

• The probability that all tracks in the event come from the primary vertex was 
required to be less than 10%, as described in [22]. The events were then considered 
as 6-fiavour candidates. This choice is the result of a compromise between the signal 
yield and the level of combinatorial background in the A+ mass window; it gives a 
6-tagging efficiency of 80% with a purity of 54%. 

• The sign of the lepton charge must be opposite to the A+ charge. It can be seen 
in Figure 1 that there is no evidence for a A+ signal in the wrong charge mass 
distribution. The upper limit measured from the data using the wrong sign lepton 
rate is 7% (95% C.L.) of the right sign sample. Similar limits have been obtained 
from a simulated sample of 66 events: a fit to the A+ mass distribution of A+Z~ 
candidates, as in the actual data analysis, gave 438 A^Z~ combinations. A sample 
of 13 ± 5 A+ candidates were associated to fake leptons. Their fraction is then 
(4.3 ± 1.5)% of the signal ^ 

All final states studied in the following must have a track satisfying the lepton selection 
requirements. 

^This number has been corrected for the different rate of fake leptons in the simulation and in the real data, and it has 
been subtracted from the amount of candidate events in the estimates of absolute rates. 
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2.1 2.2 2.3 2.4 

Mass of Ac all channels (GeV/c^) 

Figure 1: T/ie A+ mass distribution including the decay channels pK^n^ , Att+tt+tt^ and 
pKg: A^l~ combinations (white histogram) , A^/+ combinations (hatched histogram) . The 
curve is the sum of the fitted functions for each channel (see section 4-6)- 

4.3 A+ pK 7r+ selection 

Triplets of charged particles of total charge unity, each track with at least one hit in 
the microvertex detector were selected. The momenta were required to be larger than 
3 GeV/c (proton candidate), 2 GeV/c (kaon candidate) or 1 GeV/c (pion candidate), and 
the total momentum to be larger than 8 GeV/c. The mass of the A+ candidate had to 
lie in the 2.1-2.49 GeV/c^ range. A secondary vertex was fitted with these three tracks, 
requiring a probability larger than 0.001. The primary vertex was found iteratively 
using initially all the tracks of the event and following the procedure used in [26]. The 
lepton track was not a priori excluded from the general vertex fit to avoid a possible bias 
on the lepton side. The combined A° — A+ flight distance was then computed as the 
difference between the secondary and the primary vertex. It was signed with respect to 
the momentum direction of the triplet, and its projection on the plane transverse to the 
beam axis had to be larger than +0.02 cm. The particles compatible with both p and K 
identifications were kept, and no identification was applied to the pion candidate. The 
reconstruction efficiency for this channel is 8%. 

The invariant mass distribution is shown in Figure 2a. The curves in Figure 2 are 
obtained by fits of Gaussian distributions in the signal region, added to a linear back- 
ground. The Ac mass is measured in the pK~7i~^ channel, and fixed to the same value in 
the other two. The resolution a of the gaussian is left free in the pK^n^ channel, and 
found to be 13.5 ± 1.8 MeV/c^. The number of A+ pK-7i+ events is 80.4 ± 15.0. 
In the other channels, discussed in the following sections, the resolution was fixed, and 
derived from the previous one according to the ratio of the simulated resolutions. After 
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reconstruction in the simulation, the gaussian resolutions are found to be 12 MeV/c^ for 
the pK~7r~^ channel, 15 MeV/c^ for Att+tt+tt" and 16 MeV/c^ for pK^. 
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Figure 2: A^ mass spectrum in the a) pK tt^ channel; b) Att+tt+tt channel; and c) 
pKg channel. The fitted functions are described in section 4-3. 



4.4 A+ Att+tt+tt selection 

The A-hyperon candidates were selected by the DELPHI algorithm which uses the 
presence of a remote decay vertex to tag the A —> p7i~ candidates, as described in [22]. 
The hyperon and A+ momenta were required to be larger than 2.5 GeV/c and 10 GeV/c, 
respectively. All three pions were required to have a momentum larger than 0.4 GeV/c, 
and their tracks to have at least one associated hit in the microvertex detector. The charge 
of the triplet had to be +1 for A and -1 for A, and the three pions had to form a common 
vertex (the A'^ decay vertex) with a probability larger than 0.001. The projection of the 
flight distance of the A+ transverse to the beam had to be larger than +0.02 cm. 

There is some evidence for a A+ signal in the ATr+vr+vr" invariant mass distribution, 
as shown in Figure 2b, with a fitted signal of 15.1 ± 5.7 A^ Att+tt+tt" candidates. 
All candidate combinations are shown in this figure. Some events contain more than 
one accepted combination, and for two events, both of them are in the A+ mass range 
2.260 < M(A+) < 2.310 GeV/c^. 
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4.5 A+ ^ pKg selection 

Once a candidate had been found [22], an identified proton was searched for, 
such that the Ac candidate would belong to the hemisphere defined by the lepton track. 
The momenta of the and proton had to be larger than 3 GeV/c. The {pK^) mass 
distribution is shown in Figure 2c, with a fitted signal of 7.9 ± 4.0 candidates. 

4.6 The full A+l" sample 

The ratios of the number of events obtained in the three channels are compatible with 
the known A+ branching fractions [25]. The curve in Figure 1, which shows the full A+Z~ 
sample, is obtained from the sum of the three Gaussian functions, with their widths as 
indicated in section 4.3, and a unique central value. The Gaussian functions are weighted 
by the number of observed events in each channel, and with a free overall normalisation. 
The number of A+Z~ candidates is found to be 101.2 ± 17.4. 

5 Other sources of A+l~ final states 

There are several physical processes which lead to the A+-lepton final states, in addi- 
tion to the signal from A^ — > A'^l^ui, named elastic channel in the following. They arise 
from instrumental or physical backgrounds: 

a) the semileptonic decay A° A'^T~Ur with the subsequent decay r~ /"z/^z/^; 

b) fake leptons: the fraction of fake leptons, estimated to be (4.3 ± 1.5)% in section 4.2 
will be subtracted from the sample of A+Z" candidates in sections 7 and 8; 

c) decays of B mesons: B — > A+iV/~P;(X), where N is an antibaryon and X stands 
for any number of or (tt+tt") pair; 

d) decays into other charmed baryons A° A*+/^P/ , or A° (ScTt)"^/^//;; 

e) non- resonant A^l^uiX final states. The multiplicity Nc of charged particle tracks 
compatible with the A° vertex, defined as the combined (AJ) vertex (ignoring the 
Ac lifetime), the missing mass in the A° decay, and the A^/~ mass will be used to 
investigate this component, including as well the decays from d). The sum of d) and 
e) final states is named inelastic channels in the following; 

f ) other weakly decaying 6-baryons, such as Ej,. The production rate of S& from b quarks 
is however 10 times lower than the Af, production (as discussed in section 7), and 
the fraction of Ac final states from S5, if assumed to be similar to that of Sc decays 
into A baryons, is also less than 10%. This background has thus been neglected. 

The final states b) to f) are background processes, while the r lepton final state a), 
belongs to the signal. Its contribution was estimated with the full simulation of the decay 
A° —>■ A+T~z/i-, assuming the same couplings as for the other leptons. The branching 
fraction is found to be six times smaller than for A'^l~Pi with light leptons. The observed 
T contribution is further suppressed by the lepton momentum selections. The estimated 
yield amounts to 2 events in the full A^ sample of 101 events and can be neglected. 

The number of events from the decay of B mesons in the A^ sample was computed 
with the full simulation, using the value quoted in [25]: Br{B A+7VZ~i/;) < 0.0032 at 
the 90% C.L., which relies on the data from [27]. The contamination of such decays into 
the sample of A+/~z/; events is found to be less than 1.5 events (90% C.L.). The inclusive 
final states with an additional charged or neutral pion have a branching fraction smaller 
than 0.0064 [25], and their acceptance is found from the simulation to be smaller by a 



8 

factor 0.62, so that they can contribute up to 1.9 events. The total contamination of 
Iv^ N Xl~ Di events from B mesons into the A+/~ sample is thus estimated to be less than 
3.4 events (90% C.L.). 

6 Charmed baryon contributions to b-baryon decays 

To extract the form factor corresponding to the A° — > A+/~z/; decays from the A^/~ 
sample, the contributions arising from the elastic and inelastic channels must be evalu- 
ated. Whether resonant charmed hadrons are present in the corresponding mass spectra 
is investigated first: if there is a large and dominant resonant contribution, its produc- 
tion should be described by the appropriate form factors. As shown in the following, no 
signal is observed and the corresponding upper limits will be given. Other experimental 
distributions which are sensitive to the fraction of inelastic channels in the sample are 
also considered. As it will be shown, they can only be understood under the assumption 
of a substantial contribution from inelastic A^tttt final states. 

6.1 The resonant states 

In addition to the elastic A+Z~P/ channel, several charmed baryon final states can 
contribute to the A+/~ sample, such as (Ec7r)"'"XZ~P/, and A*+XZ'~P/, in which the Ec and 
A*"*" decay into a A+. 

The isospin of the hadronic final state should be / = within HQET. In this section, the 
presence of the HQET-allowed hadronic final states which correspond to the following 
decay channels, all with a A+tttt final state, are investigated: 

• AO ^ E++7r-rP,, with E++ ^ A+7r+; 

• A" SOvr+Z-P;, with E" A+vr-; 

• AO ^ E+7rO/-P/, with E+ ^ A+vr"; 

• AO ^ A*+/-Pi, with A*+ ^ A+TT+TT^, or A+ttOttO. 

The search for resonant states described in this section is not sensitive to channels with 
TfO's. The first three decay modes are expected to have the same branching fractions. 
The Ec(2455) and Ec(2520) resonances would show up as peaks in the distribution of the 
variable Qs, defined as: 

Qs = M(A+7r) - M(A+) - m,. (8) 

As the relative sign between the A+ and the charged pion has not been distinguished, 
the Qy, distribution should contain the same number of events from A^ — > E^+tt^/^P;, 
and AO — > EOtt+Z^P;, and the same combinatorial background. For display purposes, 
the expected spectrum for the decay AO — > (Ec(2455)7r)/"Pi is shown as a dashed line in 
Figure 3, under the assumption that the observed number of events satisfies: 

K,,{Kl ^ A+(X)Z-pO 

where the denominator is the total number of events in the A+Z~ sample. The Ec(2520) 
signal, indicated as a dotted line histogram, is broadened due to its natural width of 
20 MeV/c^. 

From the number Nobs of events found in the signal region (column 3 of Table 1), 
upper limits have been derived on the number of observable decays from the A^ — > 
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Figure 3: Distribution ofQ^, = M(A^7r)— M(A^)— in data (solid line histogram) . The 
simulated T,c(2455) (dashed line histogram) and Sc(2520) (dotted line histogram) signals 
each assume NoUA'l ^ S++7r-/-p0 + Ar„6,(A0 SOyr+Z-pO = 0.1 Nobs{At{X)l-iyi). 
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Figure 4: Distribution of Qa = M(A+7r7r) — M(A+) — 2mT, in data (solid line 
histogram). The simulated A* signals (dotted and dashed line histograms) assume 

KUK ^ a:+/-z70 = 0.1 N,UKix)i~^i)- 
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Decay channel 


window 

(^LTeV/c ) 


A T 

Nobs 

observed 


A T 

estimated 


limit 

upper 


1 • 'J- 

limit on 

-^channel 


^ Ee(2455)++7r-/-z/i + 
AO ^ Sc(2455)°7r+/-z>i 


0.025-0.030 


3 


2.5 


5.30 


0.078 


Aj; ^ E,(2520)++7r-/-z/, + 

A^ ^ E,(2520)'V+/-/>/ 


0.080-0.105 


15 


10.5 


12.5 


0.190 


a;; ^ Ae(2593)*+/-i/, 


0.024-0.036 


1 


0.5 


4.30 


0.064 


A^ ^ Ae(2625)*+rPi 


0.054-0.066 


1 


1.5 


3.20 


0.048 



Table 1: Qb window (B = E, h.), number of observed events Nobs, estimated number of 

background events Nbkg, 95% C.L. upper limit on the observable number of signal events 
Numit, (ind the upper limit on the channel contribution to the observed sample Rchannei- 



(EcTt)"'"/ Pi channels (Table 1). These limits are obtained from the highest expected yield 
of signal events, TVs, such that the probability: 

ProbiNj: + Nbkg < Nobs) > 5%. (10) 

The probability law for the number of events is assumed to be a Poisson law, with 
the mean given by the sum of the background contribution Nbkg, interpolated from the 
adjacent bins of the Qs spectrum, and N-^, the expected yield. The uncertainty on the 
estimated mean background level has not been included. 

To search for the A*"*" states, the distribution of the variable 

Qa = M(A+7r+7r-) - M(A+) - 2m^ (11) 

is considered (see Figure 4), and 95% C.L. upper hmits have been obtained using the 
same procedure as in the (Ett) final states. The kinematical window for Qa is given by 
the column 2 of Table 1. 

The upper limits in column 5 of Table 1 refer to the observable final states. The con- 
tributing fraction Rchannei columu 6 is: 

jj N channel ^10^ 

^ Nobs{^t{X)l-ui) ^ ' 

where Nchannei is the (maximal) number of events from this channel contributing to the 
A|f/" sample: a correction factor of 3/2 has been applied to the upper limits for the A*+ 
final states, as the pion pair is assumed to have / = 0. The same factor 3/2 has been 
applied to the observed E+ contribution in order to include the A+7r° final state. 

6.2 The A+TTTT contribution 

No evidence for the production of excited charmed states in A° semileptonic decay 
has been found, but the limits obtained still allow for a substantial contribution from 
these final states. As, in addition, non-resonant charmed channels might be present, an 
inclusive approach will be adopted to evaluate the combined contribution of the resonant 
and non-resonant final states. This inelastic contribution is investigated by studying 
three distributions: 
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• the multiplicity of charged particle tracks associated to the secondary vertex 
(Figure 5a); 

• the missing mass squared M^^^^ = {p^^ — Pj^+ — ViY (Figure 5b); 

• the A+-lepton mass, which is expected to be smaller for A^(7r7r)/~z/; final states 
(Figure 5c). 




0123 024 3 4 5 

Multiplicity Missing mass squared Mass (Ac lepton) (GeV/c^) 



(GeV/c^)^ 

Figure 5: The charged multiplicity; b) the missing mass squared and c) the (AJ) mass 
distributions in the A^ mass window. The hatched areas show the A^7r"'"7r~ contributions 
(including ). The cross hatched histogram is the combinatorial background. 



6.2.1 The multiplicity 

The multiplicity of charged particles other than the A+ decay products and the lep- 
ton, compatible with the A° vertex was evaluated by a neural network algorithm which 
separates tracks compatible with the primary or the secondary vertex. The probability of 
the vertex assignment was required to be larger than 0.5. The b and c vertices were not 
separated in this treatment. The result for the events situated in the A+ mass window 
2.260 GeV/c^ < Ma, < 2.310 GeV/c^ is shown in Figure 5a. The contribution of the 
A+7r°7r°/~z/i channel is assumed to be 1/2 of the charged pion mode A+tt'^'tt^/^z/;, as the 
hadronic state is assumed to have / = 0. 
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J-'miss ~r ■'-'same — „~r n r~ ' \ ) 



6.2.2 The missing mass in — > A+1 z/i 

To reconstruct the missing mass squared M^^^^, with M^-^^ = {p^^ — — PiY, the 
four-momentum of the was evaluated assuming the decay channel A° A'^l~i?i. The 

energy is = E\^ + Ei + Ey, and the A^ momentum has been computed as E\^ — 

with = 5.624 GeV/c^. The neutrino energy and the direction of the have to be 
determined. The energy of the undetected neutrino in the exclusive A^ — > A^l^j?i channel 
was estimated as in previous studies of semileptonic decay of B mesons [1]. The total 
energy in the A+/~ hemisphere, which is the sum of the visible energy (Esame) in the same 
hemisphere, and of the neutrino energy (Emiss), was computed using the observed total 
masses of the A+ hemisphere (Msame) and of the opposite hemisphere (Moppo), together 
with four-momentum conservation: 

V"' I same o; 

As Esame, ^samei ^nd Moppo, are approximately known due to the detector inefficien- 
cies, an empirical correction fsim{Esame) (as in [1]) has been estimated from a simulation 
of the exclusive semileptonic channel, to improve the accuracy on the neutrino energy 
reconstruction: 

Whenever this procedure leads to a negative energy, the value E^j^ = is used. The final 
resolution on the neutrino energy is around 33%. Adding this energy to the energy of the 
A+ and of the lepton gives the energy of the A^ with a resolution of 8% (Gaussian fit). 
The direction of the A° can be estimated by two different methods: the momentum of the 
A+/~ system gives the best accuracy when the line of flight is short, while the positions of 
the primary and secondary vertices, work best at large separations. In the data from 1992 
and 1993, where the microvertex detector did not provide any z information, the A+/~ 
momentum was always used to measure the polar angle 9. The weighted combination of 
the two which gave the best resolution on the Isgur-Wise variable w was chosen. 

It is seen in Figure 5b that the reconstructed missing mass squared M^-^^ is sensitive 
to the presence of the inelastic channel. 



6.2.3 The A+1 mass 

The A+/~ mass distribution is expected to be shifted to lower values in the inelastic 
A+7r7rZ~P; final states. This effect is apparent in Figure 5c and it constrains the proportion 
of inelastic events 



6.2.4 Fit of the elastic fraction fA^ 

To measure the fraction /a^ of elastic A+Z~P/ decays, the three previous distributions 
are considered to be the sum of three components: 

• the A° A+/~z/; decays; 

• the A^ — > A+7r7rZ~P/ decays. As this channel is a sum of many states, it is simulated 
without a form factor, using the quark matrix elements for weak decays in PYTHIA 
[28]; 

• events from the combinatorial background present under the A+ signal. The shape 
of this component is evaluated using events situated in the side bands of the Ac mass 
peak. Its normalisation is fixed according to the fit of Figure 1. 



13 



The overall normalisation is fixed to the total number of observed events. The elastic 
fraction /a^ in the final state is obtained from three fits to the three distributions in 
Figures 5. Each of them is first adjusted independently. The statistical correlations of 
the three results are obtained from the simulated two dimensional distributions of the 
three pairs of variables. They are respectively: 0.25 for (multiplicity, A^/^ mass), 0.35 
for (multiplicity, missing mass), and 0.65 for (missing mass, A+/~ mass). The combined 
fit to the three values gives: 

The ratio /a^ is not a branching fraction, but the ratio of the observed contributions in 
the data sample. The systematic uncertainties arising from the identification efficiency, 
the time of flight and the modelling of the 6-quark fragmentation are negligible. The 
uncertainty on fP, as measured in this paper, changes by ±0.02, and the uncertainty 
from the simulation of the inelastic channels is estimated to be ±0.03 by substituting 
a EcTT decay for the A+tttt non-resonant prediction. The systematic uncertainty from 
the combinatorial background is estimated by comparing side bands of different sizes and 
contributes ±0.05. The limits on the observed resonance contributions given in section 6.1 
can be turned into limits on branching fractions. The fraction of observable decays can be 
derived from the assumption of an / = final state and is found, as already mentioned, to 
be 2/3. The kinematical acceptances for (Scvr, A*)/~t/; and A+/~t/i in the charged decay 
channels differ, with e(A+(7r+7r~)Z~P;)/e(A+Z~P;) = 0.57, according to the simulation. 
The branching fraction to the exclusive A.'^l'Pi final state can then be estimated to be: 

Although half of the A+/~ events in the A° decay arise from a K^'K'kI'vi final state, no 
signal from resonant production of Ec or A* has been observed in the present sample. 



7 Leptonic branching fraction from the event rate 

The previous determination of the fraction of clastic scmilcptonic decays of A° 
K^l~vi in the A+/^ sample allows a measurement of the semileptonic branching fraction, 
Sa^, in this transition. As the total width of the A° is given by its lifetime, this branching 
fraction can then provide a direct measurement of the slope parameter of the form 
factor. The branching fraction i?Ac has been measured from the number of A° — > A+Z~P; 
candidates using only the decay channel A+ — > pK~T:^ , which has less background: 

N{Ktri>i) ^N^2R, f{b ^ A°) B^^ Sr(A+ ^ pK-7r+) eipKnl) Ni. (17) 

As the event rate is much more sensitive than the w shape to the knowledge of the 
identification efficiencies, this analysis has been restricted to the periods of data taking 
where both gas and liquid RICH were present. The number of hadronic Z decays is 
iV| = 1.52 X 10^ and the fraction of ?^flavoured final states is Rb = 0.217 [25]. Detailed 
comparisons have shown that the identification efficiencies in data and simulation are 
then in excellent agreement. The simulated reconstruction and identification efficiency 
is found to be, including the Icpton, e{pK7rl) = (8.0 ± 0.8)%. The branching fraction 
5r(A+ -^pK-n+) = (5.0±1.3)%, and /(6 ^ A^) = 0.108±0.020 (subtracting 1% for Eb 
from the quoted value in [25]) are used. The number of lepton families is Ni = 2, as the 
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rererence value 


± uncertainty range 


SB^JBa^ 


Br{Aj pK 7r+) = 0.050 


0.013 


0.26 


/(6 ^ Afe) = 0.108 


O.OzU 


O.lo 


/a. = 0.60 


0.10 


0.16 


e(pir7r/) = 0.080 


0.008 


0.11 


ta, = 1.23 ps 


0.08 ps 


0.065 


B meson decays < 3.4% 


0.034 


0.034 


fake leptons = 4.3% 


1.5% 


0.016 


Total {6BaJB^J 




0.37 



Table 2: Main sources of systematic uncertainties on the expected rate. 



T contribution is negligible. The total number of observed A+/~z/; and A^Xl^Ui events 
(without the enrichment selections for the elastic channel of section 8) is 47 ± 10. The 
fake lepton contamination, as determined in section 4.2 amounts to 2.0 ± 0.7 events and 
should be subtracted. The "observed" number iV(A+/"z/i) of exclusive A+Z~P/ decays can 
be estimated from the latter fraction to be A^obs = /ac " (47 — 2.0) = 27 ± 6. This imphes 
i?Ac = (4.7 ± 1.1 (stat))%. The main sources of systematic errors on B\^ are given in 
Table 2 and sum to a relative error of 37%, so that: 

Ba, = (4.7 ± 1.1 (stat) ± 1.7 (syst))%. (18) 

Most of the systematic error arises from the uncertainty on -Br(A+ — > pK~n'^). The other 
systematic errors include the uncertainty on the reconstruction efficiency e, including 
identification, the contribution of the 5-meson decays which was varied up to the maximal 
value of 3.4% found in section 5. The impact of the lifetime uncertainty on the acceptance 
is negligible, but its effect on the parameter via the normalisation of the branching 
fraction is large and it is included here to simplify the presentation. The contribution of 
the pK~Tr~^ branching ratio can be explicitly extracted: 

= (4.7 ± 1.1 (stat) ± 1.3 (Br(Ac pK-7r+)) ± 1.3 (other syst))%. (19) 

The semileptonic decay width can be computed from p^, using equation (2), under 
the assumption that ^b(I) = 1. The total width is given by the A° lifetime ta^ = 
1.23 ±0.08 ps [25], and the semileptonic branching fraction B\^ provides (within HQET) 
an estimate of the slope fP: 

pI,,, = 2.05l°i° (stat error only) . (20) 

This measurement of p^ will be combined in the next section with a fit to the distri- 
bution of the Isgur-Wise variable w to obtain an improved determination of the slope 
parameter and of the branching fraction. The impact of systematic errors will be evalu- 
ated for this combined fit. 

8 Combined fit to the w shape and the event rate 

8.1 The enriched A+-lepton sample 

As has been shown in Figure 1, the number of A^-lepton pairs is 101 ± 17, and this 
sample contains the elastic A^l~ui and inelastic A^{mT)l~ui final states. The charged 
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multiplicity N^., as well as the (A+ lepton) mass can be used to enrich the sample with 
respect to the elastic channel, by selecting multiplicities smaller than 2, and A+/~ masses 
larger than 3.5 GeV/c^. The number of A+ candidates left is obtained from a fit to 
the mass distribution of the candidates, shown in Figure 6, and is 62.5 ± 10.0. The 
remaining AcTttt contribution is obtained from the simulated efficiency, after application 
of the enrichment selection, and is found to be 10.0 ± 2.8 events. The fraction of elastic 
A.^l'Pi events in this sample is r^^ — 0.84±0.17. The Ac mass distribution of the enriched 
sample, together with the different components is shown in Figure 6. The slope parameter 
of the w distribution is determined from this enriched sample. The candidates obtained 
from the three decay modes of the A+ considered in section 4 have been separately 
analysed. In each of these samples, the non-A+ background is measured from a fit to the 
A+ mass distribution. 

The complementary sample, enriched in A+tttt, will be used to monitor the w-shape 
of the A+TTTT background. 

8.2 The w-shape likelihood 

The four-momentum of the meson is reconstructed as described in section 6.2.2. 
The values of and w can then be estimated. The resolutions achieved are similar to 
those obtained in [4], and Aw/w is close to ±8%. 

A likehhood fit to the two-dimensional (M(Ac), w) distribution is then performed, with 
2.190 < M(Ac) < 2.385 GeV/c^ and 1.0 < w < 1.6. The mass and w dependences are 
assumed to factorise in the probability distribution for each event k in channel i, where 
the index i = 1,2,3 runs over the three input channels (three final states, the two lepton 
samples are combined): 



The Gaussian term G{Mk) describes the A+ contribution, "{w) is the inelastic 
background, while the combinatorial background B^"^" (w) is mass independent in the Ac 
mass window investigated, as justified by a direct inspection of its shape in the simulation 
in Figure 6. The coefficients // (for the elastic Ac signal), f^"'^^, and ff are fixed, and 
are obtained from the pKir mass spectrum and the elastic fraction of 0.84. 

In each channel i, the contribution of the signal Si{w) is obtained from the simulated 
events weighted at a given by the squared form factor: cxp{—2p^{wy — 1)) where Wg 
is the generated value of w. The function Si{w,p'^) is a convolution of the physical w 
distribution with the detection efficiency and the resolution of the reconstruction of this 
variable. It has been expressed as a factorised expression, as it was easier to parametrise 
directly Si rather than the resolution function: 



The function Si{w,p'^ = 0) is parametrised as Si{w,p'^ = 0) = a{w — l)™'exp(p(w — 1) + 
q{w - 1)^)), with a = 0.458, m = 0.154, p = 17.2, q = 40.3. A very good description of 
the simulated w distributions was obtained when the exponent a{w) was assumed to be 
a linear function of w: 



Si{w,p'^) = Si{w,p'^ = 0) exp(p^ • a{w)). 



(22) 



a{w) — ao + CKi • {w — 1). 



(23) 



The coefficient ao is a normalisation coefficient, while ai describes the w dependence of 
the signal. In the absence of smearing and detector effects, ai would be -2, given the 
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Real data 

bckg from MC 

□ AcX 

□ bckg from fit 




C-i I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

2.2 2.24 2.28 2.32 2.36 

Mass of Ac(GeV/c2) 



Figure 6: The mass distribution in the (enriched/elastic) sample o/A+/~ events used in 
the likelihood shape analysis (all channels). The estimated background (dotted histogram) 
is obtained from simulated hadronic Z decays and normalised to the actual data 



exponential parametrisation in equation (7), and the actual value ai = —1.67 ± 0.08 is 
close to this estimate. 

The Bf^iw) and B^°^''{w) functions are found from the data. The shape of B^^iw) is 
obtained from the w distribution of the subsample of events enriched in A+tttt, with 
a charged multiplicity at the A° vertex N^. > 2 oi with a (A^/~) mass < 3.5 GeV/c^. 
A correction factor derived from the simulation is applied to the w distribution of this 
sample, to relate the background shapes in the enriched (elastic) and anti-enriched (AcVrvr) 
samples. This background, corrected for the enrichment bias is then described by a 
function B^^{w) = {w-l^^ exp(-6i(w - 1)), with ai = 2.22 ±0.50, and hi = 27.5 ±5.8. 
The non-A+ background, B^"^''{w), is evaluated from side-bands in the mass spectrum, 
chosen outside the mass window 2.260 < M(Ac) < 2.310 GeV/c^. Its shape is described 
by the same parametrisation as the previous one, with new values of a and b. In the 
dominant pK~7r~^ channel, for instance: a = 4.28 ± 0.34, and b = 22.3 ± 1.9. 

The distributions of the signal in the different channels Si{w,fP), and the shapes of 
the backgrounds Bj{w) are normalised to unity, and the coefficients ff , have been 
measured from the data and are fixed. The likelihood is the sum: 

jC-shape = — J2k Log(P(Mfc, Wk)), and the slope obtained from the one parameter fit to the 
w distribution is: 

= 1.59 ± 1.10 (stat). (24) 

The quality of this fit can be checked in Figure 7a, where the predicted distribution has 
been normalised to the observed number of events in real data. The /NDF between 
the distribution predicted from the likelihood fit and the observed distributions is 5.2/11. 
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8.3 The combined event rate and w shape hkeUhood 

The information on the shape and on the absolute rate can be combined into an 
optimised determination of p^, assuming ^b(I) = 1, and a new hkehhood fit is performed 
where the observed number of events is included as a constraint. The expected number of 
events N^.^ is derived from the semileptonic branching fraction, which is itself a function 
of as described in section 7: 

C- = Cshape - Nobs'i^Og{Nex) + N^x- (25) 

The full statistics are used for the shape likelihood, while as in section 7, the rate Nex is 
measured only when the gas and liquid RICH are simultaneously operational, and for the 
A+ — > pK^TT'^ channel alone ^. The actual likelihood function includes the contributions 
from the different backgrounds, and the value of obtained from this one parameter fit, 
shown in Figure 7b, is: 

p2 = 2.03 ±0.46 (stat). (26) 

The statistics in Figure 7b are smaller than in 7a as only data with an operational RICH 
are used. 




Measured w Measured w 

Figure 7: Comparison between fitted and measured w distributions inside the mass 
window: a) fit to the w distribution; b) combined fit to the w distribution and the A^l~iJi 
rate (restricted to the working RICH sample). The statistics in Figure 7b are smaller 
than in 7a. 



The semileptonic branching fraction corresponding to this value of p^ is: 

Br{Al ^ Ap-ui) = {5.0tli (stat))%. (27) 

The expected number of events is 27^*^4, while the observed number of elastic Ac events 
quoted in section 7, and included in the likelihood fit, is 27 ± 6. 

There are three groups of systematic errors: the errors associated to the prediction of 
the w shape (the identification efficiency, the parametrisation of the shape as a function 

"^The statistics of observed and expected events were scaled to reproduce the actual statistical uncertainty on the number 
of events. 
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source of error 


± uncertainty range 


contribution to 




0.15 


0.08 


ow/w = U.Uo 


O.OOo 


O.Oo 


Oil — — l.Ol 


U.Uo 


n no 
u.uz 


Fragmentation (< xa^ >) 
a, ai and b, bi (shape) 


0.008 


0.05 


stat. errors 


0.01 


Enriched fraction r^^ = 0.84 


0.17 


0.10 


Anti-enriched background shape 


stat errors 


0.02 


SBaJBa^ (from Table 2) 


0.37 


+U.7U 
-1.00 


Total 




+0.72 

-1.00 



Table 3: Main sources of systematic uncertainties on the slope . 



of p^, the fragmentation), the uncertainties related to the expected yield (absohite effi- 
ciencies and branching fractions) and the systematic effects arising from the subtraction 
of the Ac background processes discussed in section 8.2 (parameters a, Oi and 6, 6i). 

The momentum dependent uncertainties are dominated by the contribution to the 
global identification efficiency of the proton and kaon identification efficiencies e^d. The 
observed momentum distributions and therefore the w spectrum are sensitive to the 
detailed simulation of the identification algorithms. The difference between the efficiencies 
of proton identiffcation in data and in simulation was monitored using a sample of selected 
A events. For each year, the identiffcation efficiency as a function of momentum was 
measured for data and simulation. The difference typically reaches 15%. Simulated events 
were reweighted by the ratio of the data and simulation efficiency, and the systematic error 
on fP' measured from the shift in fp' between the original and reweighted w distributions. 

The resolution on w is found to change by approximately 10% between different ver- 
sions of the energy reconstruction algorithms. The impact of this uncertainty was evalu- 
ated by degrading the resolution by 10% in the simulation. Given the small contribution 
of the A^ semileptonic decays with a r lepton, no systematic error was assigned for this 
component of the signal. 

The uncertainty on the number of Ec and A* ffnal states is included in the uncertainty 
on the fraction ta^,, equal to the statistical uncertainty on its value. The systematic 
uncertainty on the shape of the non-Ag background is evaluated by changing the shape of 
this background according to the statistical error on its parametrisation. This background 
is evaluated from the side-bands of the Ac mass spectrum. 

The main systematic uncertainties on are summarised in Table 3, where the domi- 
nant contribution of the normalisation error (as estimated in Table 2) is singled out. The 
value of fP is found to be: 

= 2.03 ± 0.46 (stat) t^^^ {Br{Kc pR-n-^)) 1°;™ (other syst). (28) 

The asymmetry in the errors on pP arises from the strong non-linearity of the relation 
between the normalisation and the slope over the large range of SB\^. The contribution of 
the branching fraction into pKtt to the uncertainty on the normalisation is given in Table 
2, and the associated systematic error on has been evaluated separately in equation 
(28). The corresponding elastic semileptonic branching fraction is: 

Sr(AO ^ A+rP,) = {b.Otli (stat)+;:^ (syst))%. (29) 
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To obtain p^, the prediction of HQET ^(1) = 1 has been assumed, and no theoretical 
uncertainty has been included in the systematic error to account for this hypothesis. A 
10% change of ^(1) would amount to a 20% change in the rate, and to 0.3 in p^. 



9 Conclusions 

A first measurement of the form factor of the A^ beauty baryon has been achieved in 
the A^ — > A+/~IY decay channel. Assuming an exponential behaviour of the Isgur-Wise 
function: 



-p^{w-l) 



the slope parameter was determined from a w-shapc analysis to be: fP' = 1.59 ± 
1.10 (stat). If the validity of the HQET relation ,^(1) = 1 is assumed, and the event rate 
taken into account, an improved determination of the slope can be obtained: 

= 2.03 ± 0.46 (stat) toio (Br (Ac pX~7r+)) to.lo (other syst). 

The evaluation of the systematic errors takes into account the actual variation of 
arising from each source. 

The semilcptonic branching fraction into the exclusive semileptonic mode was mea- 
sured within the hypothesis ^b(I) = 1 to be: 

Br{Al Atl-Pi) = (5.0l^i (statj^J;? (syst))%. 

The fraction of elastic A+Z~P; events is found to be: 

r(AO ^ Atl-U,) _ Q_^^^,10 (3,,,)+0.07 _ 



r(AO ^ A+/-Z/0 + r(AO ^ A+TiTTl-ui] 

The spectra of multiplicity, missing mass, and A+/~ mass shown in this paper strongly 
hint at such a sizeable fraction of non-elastic hadronic modes. This inelastic contribution 
is larger than assumed in the experimental determinations of the A^ lifetime, such as [29] , 
and will affect its value, as the ratio of the Ac and A{, momenta Pa^/pa^ is 20% lower in 
the inelastic channel. 

In spite of the evidence for a large inelastic contribution in semileptonic decays, no 
indication of charmed baryon resonances has been found in the final state. 

The parameter refiects the structure of the A° baryon. Its value is somewhat larger 
than in the 6-meson channel, where p^ ~ (0.60 to 1.3) as measured in [2-4]. A recent 
result on the B meson decays from [5] suggests an even higher value of the slope of the 
Isgur-Wise function for mesons, with = 1.61 ± 0.09 (stat) ± 0.21 (syst). In all models 
proposed so far [16,18,30,31], the value of p^ is expected to be larger for baryons. 
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